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Abstract 

We employ a high-resolution Köppen climate classification dataset to examine 
shifts in Tundra zones within the Alps and Asia. Our analysis shows substan- 
tial reductions in Tundra areas by the mid-21st century under different Shared. 
Socioeconomic pathways (SSP1-2.6, SSP3-7.0, SSP5-8.5). Tundra zones in the 
Alps and the Tibetan Plateau are crucial for their unique climates and role as 
water reservoirs. Characterized by short, mild summers and long, severe win- 
ters, these zones are vital for the glaciers and perennial snow. The projected 
climate instability may significantly reduce alpine snow cover by mid-century 
with irreversible consequences. A 2°C temperature increase from the 1981- 
2010 baseline could eliminate the Tundra climate in the Alps and reduce it by 
over 70% in Asia. This is particularly concerning given that rivers from the 
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1 | INTRODUCTION 


Policy-makers require high-resolution climate datasets 
for creating effective local adaptation strategies (Flato 
et al., 2014). Detailed climate data strengthen cross- 
sectoral impact assessments by enabling researchers to 
consider the complex interactions between climate vari- 
ables and socioeconomic factors (Dunford, Smith, 
et al., 2015; Jäger et al., 2015). For instance, in agricul- 
ture, high-resolution data support models that forecast 
changes in crop viability and water resource availability 
under different climate scenarios (Brown et al., 2015; 
Dunford, Smith, et al., 2015; Holman & Cojocaru, 2011). 
Climate change indicators, such as the Köppen climate 
classification (Geiger et al., 2012; Peel et al., 2007), are 
essential tools for climate impact studies, especially in 
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regions with complex topography (Brown et al., 2015; 
Dunford, Harrison, & Rounsevell, 2015; Fallah et al., 
2023). The Köppen system categorizes the climate zones 
based on temperature and precipitation. 

Applying the high-resolution Köppen system helps to 
identify slight but meaningful shifts in climate zones 
(Beck et al., 2018). The Tundra zone is characterized by a 
cold environment and short growing seasons, which are 
vulnerable to warming. The Tundra is an essential car- 
bon reservoir. Warming temperatures threaten to trans- 
form these areas into carbon sources and accelerate 
climate change. 

In this study, we use the 1-km Köppen dataset from 
the Climatologies at High Resolution for the Earth's Land 
Surface Areas (CHELSA) (Karger et al., 2017, 2020, 
2021). These high-resolution data are essential for impact 
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modeling, providing an ecologically meaningful classifi- 
cation of climate zones, which is vital for assessing biodi- 
versity responses to climate change (Brown et al., 2015; 
Holman & Cojocaru, 2011; Kebede et al., 2015). 

These data provide an ecologically meaningful classifi- 
cation of climate zones for assessing biodiversity changes 
due to climate change (Brown et al., 2015; Holman & 
Cojocaru, 2011; Kebede et al., 2015). Our analysis focuses 
on the Alps and the Tibetan Plateau, particularly empha- 
sizing the Tundra climate zone. The Tibetan Plateau, 
which is warming faster than the global average, is a criti- 
cal water source for major rivers such as the Yangtze, Yel- 
low, Mekong, and Indus (Clima & Te, 2023). The Tibetan 
Plateau is among the essential drivers of the atmospheric 
dynamics in the Northern Hemisphere, affecting monsoon 
systems and weather patterns across Asia (Fallah 
et al., 2016; Song et al., 2010). Melting glaciers jeopardize 
water supplies and increase the risk of floods. As perma- 
frost thaws, greenhouse gases are released, further intensi- 
fying global warming and destabilizing landscapes. 

The Alps are critical for rivers like the Rhine, Rhone, 
and Po. Understanding how Tundra zones are likely to 
shift allows for better management of natural resources, 
protection of ecosystems, and preparation of communi- 
ties for future scenarios. 

We address the following scientific questions: 


e Q1: Which climates might replace the Tundra climate 
zones in the Alps and Asia in the future? 

e Q2: What are the projected impacts of surpassing a 2°C 
temperature increase on the geographical distribution 
of the Tundra climate zone? 


2 | DATA AND METHODS 


2.1 | Shared socioeconomic pathways 
The Shared Socioeconomic Pathways (SSPs) are scenarios 
developed to study the potential impacts of climate 
change under different socioeconomic conditions. These 
scenarios are utilized with Representative Concentration 
Pathways (RCPs), which describe different trajectories for 
greenhouse gas concentrations. The Intergovernmental 
Panel on Climate Change'’s (IPCC) Sixth Assessment 
Report (AR6) extensively reviews literature using SSP 
scenarios outcomes (Eyring et al, 2016; O'neill 
et al., 2016). 

Here, we present a summary of the SSPs used in this 
study: 


e SSP1-2.6: Represents a sustainable development sce- 
nario that maintains global warming below 2°C. This 
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TABLE1 Details of the CMIP6 models used in the study. 
Model Resolution Components 
GFDL-ESM4 Atmosphere: Atmosphere: FV3 

1° x 1° Ocean: MOM6 

Ocean: 1° x 1° Land: LM4.1 
IPSL-CM6A- Atmosphere: Atmosphere: LMDZ6 
LR De xe D Ocean: NEMO v3.6 

Ocean: 1° x 1° Land: ORCHIDEE 
MPI- Atmosphere: Atmosphere: 
ESM1-2-HR 1° x 1° ECHAM6.3 

Ocean: 0.4° x 0.4° Ocean: MPIOM 

Land: JSBACH 

MRI-ESM2-0 Atmosphere: Atmosphere: GSM4 

EE Ocean: COCO4.9 


Ocean: 1° x 1° Land: MRI land model 


UKESM1-0-LL Atmosphere: 
1.25° x 1.875° 
Ocean: 1° x 1° 


Atmosphere: 
HadGEM3-GC3.1 
Ocean: MEDUSA 
Land: JULES 


scenario requires strong climate policies and mitiga- 
tion of greenhouse gases. 

e SSP3-7.0: Depicts a world characterized by regional 
conflicts, slow economic growth, and moderate to high 
greenhouse gas emissions. 

e SSP5-8.5: Represents a scenario with rapid economic 
growth dependent on fossil energy, leading to maxi- 
mum greenhouse gas emissions and significant global 
warming. 


We use the datasets base on five Coupled Model Inter- 
comparison Project Phase 6 (CMIP6) models, identified 
as “priority” within the Inter-Sectoral Impact Model 
Intercomparison Project phase 3b (ISIMIP3b): GFDL- 
ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, 
and UKESM1-0-LL (Lange, 2019). Three models have 
low climate sensitivity, and two have high sensitivity, 
ensuring a balanced representation of the CMIP6 ensem- 
ble. Furthermore, these models were selected due to their 
structural independence in ocean and atmosphere com- 
ponents. Their process representations were considered 
reliable by climate impact modelers (Lange, 2021). Addi- 
tionally, their selection aligns with previous ISIMIP 
phases, maintaining consistency and comparability with 
previous studies. Table 1 presents the descriptions of the 
CMIP6 models used in the study. 

The ISIMIP3b bias adjustment process utilizes statisti- 
cal downscaling and bias adjustment techniques to make 
model outputs like the observed ones. The bias adjust- 
ment and statistical downscaling codes of ISIMIP3BASD 
and W5E5 (Water and Global Change Forcing Data for 
the 5th version, v2.0) are employed to produce climate 
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projection input data for CHELSA (Lange et al., 2021). 
We consider uncertainty using the range of projections 
produced by different models and scenarios (SSP1-2.6, 
SSP2-4.5, and SSP5-8.5). 


2.2 | CHELSA climatologies for 
historical period 


The foundational data for CHELSA climatologies were 
derived from the European Centre for Medium-Range 
Weather Forecasts (ECMWF) Reanalysis Interim 
(ERA-Interim) reanalysis (Dee et al., 2011), which pro- 
vides large-scale atmospheric variables at a coarse resolu- 
tion of approximately 0.75° latitude/longitude. The 
downscaling process involved statistically adjusting these 
coarse-scale outputs to a finer spatial resolution of 30 arc 
seconds (~1 km?). This enables impact modelers to study 
the small-scale topographic and climatic variations essen- 
tial for ecological and environmental applications. 

The downscaling of temperatures utilized the vertical 
temperature gradients (lapse rates) from the ERA-Interim 
data, which were corrected for elevation using a digital 
elevation model (DEM). The lapse rates were interpo- 
lated to sea level and re-projected onto the terrain using 
the DEM. This allows for accurate temperature represen- 
tation across complex terrains (Karger et al., 2017). 

For precipitation, the downscaling incorporated oro- 
graphic predictors such as wind fields, valley exposition, 
and boundary layer height to refine the precipitation dis- 
tribution. The resulting precipitation fields were further 
bias-corrected using gauge-derived products from the 
Global Precipitation Climatology Centre (GPCC) 
(Schneider et al., 2014) and the Global Historical Clima- 
tology Network (GHCN) (Lawrimore et al., 2011). 


2.3 | CHELSA climatologies for future 
scenarios 


The CHELSA climatologies for future scenarios apply 
similar downscaling and bias correction methodologies 
used for the historical climatologies but are specifically 
adapted to future climate model outputs. The raw outputs 
from the CMIP6 models, which typically have coarser 
resolutions (around 1° or more), are downscaled to a 
finer resolution (30 arc seconds) using statistical down- 
scaling methods similar to those used for historical data. 
As with historical data, bias correction is applied to the 
downscaled CMIP6 projections. This process involves 
comparing the model outputs to observed data, such as 
those from the GHCN, or high-resolution reanalysis data 
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like ECMWF Reanalysis 5th Generation (ERA5). This 
comparison is used to correct systematic biases in the cli- 
mate model outputs. 

After downscaling and bias correction, monthly cli- 
matologies for future periods (e.g., 2021-2040, 2041- 
2060, and 2061-2100) are generated. These climatologies 
provide high-resolution projections of future climate 
under different SSP scenarios. 


2.4 | Validation of CHELSA 
climatologies 


Karger et al. (2017) show that validating historical clima- 
tologies involves multiple steps. First, a cross-validation 
approach is employed on the meteorological station data, 
where 20% of stations are randomly omitted, and the 
remaining 80% are used for bias correction. This process 
is repeated multiple times, showing that the bias correc- 
tion significantly improves the accuracy of the precipita- 
tion and temperature estimates. The final climatological 
products are compared against other gridded products 
and independent meteorological station data. They con- 
clude that CHELSA's predictions are comparable to or 
even better than other high-resolution climate datasets, 
particularly in complex terrains. 

The models used in these scenarios are often vali- 
dated against historical climate data using a process 
known as hindcasting. Historical inputs drive the models, 
and their outputs are compared with observations. This 
comparison helps identify and correct biases in the 
models. The iterative refinement of these models ensures 
that they better capture the vital climatic processes and 
increases the reliability of their future projections. How- 
ever, it is essential to acknowledge that uncertainties 
remain in future climate projections, and these uncer- 
tainties must be carefully considered. 


2.5 | Köppen-Geiger climate 
classification 


The Köppen-Geiger climate classification dataset was 
constructed by applying the Köppen-Geiger classification 
algorithm to high-resolution temperature and precipita- 
tion data derived from the CHELSA climatologies, 
thereby categorizing global climates based on these key 
climatic variables. 

For our analysis, we utilized the ‘kg02’ variable from 
the CHELSA bioclim dataset, a high-resolution set of cli- 
matological information tailored for ecological applica- 
tions. This layer represents the Köppen-Geiger climate 
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Köppen Climate Classification Mapping. The variables used are as follows: Tm: Mean monthly temperature (in degrees 


Celsius). Ta: Mean annual temperature (in degrees Celsius). Pa: Annual precipitation (in millimeters). Ps: Precipitation of the wettest month 


(in millimeters). P,,: Precipitation of the driest month (in millimeters). P4: Precipitation of the driest month in summer (in millimeters). 


Thmax: Temperature of the hottest month (in degrees Celsius). Tmax: Temperature of the warmest month (in degrees Celsius). T mmin: 


Temperature of the coldest month (in degrees Celsius). P;: Threshold precipitation, calculated as P; = 2 - Ta + 14 for coastal areas and 


P,=2-T, for continental areas. 


Code 
1 


AeA UOU N 


Climate classification 

Af—Equatorial fully humid 

Am— Equatorial monsoonal 

As—Equatorial summer dry 

Aw—Equatorial winter dry 

BWk—Cold desert 

BWh—Hot desert 

BSk—Cold steppe 

BSh—Hot steppe 

Cfa—Warm temperate fully humid hot summer 
Cfb—Warm temperate fully humid warm summer 
Cfc—Warm temperate fully humid cool summer 
Csa—Warm temperate summer dry hot summer 
Csb—Warm temperate summer dry warm summer 
Csc—Warm temperate summer dry cool summer 
Cwa—Warm temperate winter dry hot summer 
Cwb—Warm temperate winter dry warm summer 
Cwc— Warm temperate winter dry cool summer 
Dfa—Snow fully humid hot summer 

Dfb—Snow fully humid warm summer 
Dfc—Snow fully humid cool summer 
Dfd—Snow fully humid extremely continental 
Dsa—Snow summer dry hot summer 
Dsb—Snow summer dry warm summer 
Dsc—Snow summer dry cool summer 
Dsd—Snow summer dry extremely continental 
Dwa—Snow winter dry hot summer 
Dwb—Snow winter dry warm summer 
Dwc—Snow winter dry cool summer 
Dwd—Snow winter dry extremely continental 
ET—Polar Tundra 

EF—Polar frost 


classification and draws on the updated global climate 
zones map by Peel et al. (2007). The explanation of the 
matrices used is listed in Table 2, and the variables 
applied are as follows: 


* Tm: Mean monthly temperature (in degrees Celsius). 


EES 


Metrics 

Tm 218 C,P, > 60mm 

Tm 2 18 C,1000 > Pa > 60mm 

Tm =18 C,P, < Py 

Tm 218 C,Py <P; 

Ty <18 C,Py <10-P; 

Ta >18 C,P, <10-P,; 

Ta <18 C,Pa >10-P; 

Ta 218 C,Pa>10-P, 

—3 C<Tm<18 C,Thmax > 22 C, Ps > 30mm 

—3 C<Tm<18 C,10 C< Thmax <22 C,P;>30mm 
—3 C<Tm<18 C,Thmax <10 C, Ps > 30mm 

—3 C<Tm<18 C,Thmax > 22 C, Pa < 30 mm 

—3 C< Tm <18 C,10°C < Thmax < 22 C, P4 < 30 mm 
—3 C< Tm < 18C, Thmax < 10° C, Pa < 30 mm 

—3 C<Tm<18 C,Thmax > 22 C,Pw < P, 

—3 C<Tm <18 C,10 C <Thmax <22 C,Py <P; 
—3 C< Tm <18 C, Thmax < 10°C, Pw < Ps 

Tm < —3 C p > 22 C,P;>30mm 

Tm < —3 C,10 C <Thmax < 22 C, P; > 30mm 

Tm < — 3 C,Thmax <10 C, P, > 30 mm 

Tm < —3 C, Thmax < 10°C, P, 2 30MM, T mmin < — 38 C 
Tm < — 3 C, Thmax > 22 C, Pa < 30 mm 

Tm < —3 C,10 C < Thmax < 22 C, Pa < 30 mm 

Tm < — 3 C, Thmax < 10 C, Pa < 30 mm 

Tm < — 3 C, Thmax < 10 C, Pa < 30mm, Tmmin < — 38 C 
ENS C Ta 222 E PE 

Tm < —3 C,10 C<Thmax < 22 C,Py <P, 

TS =I Ch n < O CP SIP: 

Tm < —3 C, Thmax < 10 C, Pw < Ps, Tmmin < — 38 C 
Tere < 10) Cy nN 

Tmmax <0 C 


: Mean annual temperature (in degrees Celsius). 
: Annual precipitation (in millimeters). 
: Precipitation of the wettest month (in millimeters). 


Py: Precipitation of the driest month (in millimeters). 


Pa: 


Precipitation of the driest month in summer 


(in millimeters). 
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FIGURE1 Warming trend 
(near-surface air temperature 
anomalies with respect to 1981- 
2010) from ISIMIP3b “priority” 
models for (a) global, (b) Asia 
and (c) Alps. The Asia and Alps 
domains are shown in Figure 2. 
Boxplots show the data for 
2071-2100 vs. 1981-2010. 
SSP1-2.6, SSP3-7.0, and SSP5-8.5 
are shown in green, brown, and 
pink colors. 


* Thmax: Temperature of the hottest month (in degrees © Pr: 


Celsius). 


* Tmax: Temperature of the warmest month (in degrees 


Celsius). 


* Tmmin: Temperature of the coldest month (in degrees 


Celsius). 
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(a) Global 


IY ssps-8.5 
MW ssp3-7.0 
MN) ssP1-2.6 


Temperature Change [K] 


T T T T T T —2 
2000 2020 2040 2060 2080 2071-2100 vs. 1981-2010 


2 —2 
2000 2020 2040 2060 2080 2071-2100 vs. 1981-2010 
(c) Alps 
10 10 


—2 
2071-2100 vs. 1981-2010 


2000 2020 2040 2060 2080 


Threshold precipitation, calculated as 
P,;=2-T,+14 for coastal areas and P;=2-T, for con- 
tinental areas. 


We count the grid points labeled as Tundra and assess 
the potential shifts in the Köppen Tundra zone. 


asuaor] SUOWIWIOD aANRaID aquordde ay) Áq poussao’d ose SAIN WO ÉN jo sAN 107 Aresqry AUO KAIM UO (SuONIpUos-puL-suLIAy/WOD KoIM-AreIqrautpuo//:sdyyy) SUONIPUOD pur sua L I 99g *[PZ0T/60/0E] UO Aresqr] UUO AIIM * ZAO wWepsiog uNUEZ-z]OYWYay - tueysoy pnosep Aq 1000L2W/Z001 OI /Op/wios Aaya Aresqourpuo'sjauy//:sdyy Woy paperopumog ‘S “PZT ‘080869t1 


(b) SSP1-2.6, 2011-2040 


Meteorological Applications 


Science and Technology for Weather and Climate 


(a) Historical (1981-2010) 
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FIGURE2 Köppen climate classification for the ET (Tundra) climate zone in Asia for a) historical period of 1980-2010 and future 
scenarios of SSP1-2.6, SSP3-2.6, and SSP5-8.5 (left, middle, and right columns) and for near, middle, and far future (top, middle, and bottom 
rows). The subfigure in the top left indicates the number of grid points (each representing a 1 k m? area) that each model projected for the 
ET classification. This provides an idea of the extent of agreement among the different models regarding the size of the ET climate zone in 
the future scenario compared with the historical record. Grid points where all five models agree are only shown in the maps. 
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3 | RESULTS 


3.1 | ISIMIP3b temperature tends 

Figure 1 shows a positive warming trend over the next 
century from the ISIMIP3b climate input dataset on a 
global scale, as well as over the Asia and the Alps. The 
boxplots on the right side of Figure 1 illustrate the tem- 
perature increase during 2071-2100 compared with the 
baseline period of 1981-2010. 

SSP1-2.6 exhibits the smallest temperature increase 
across all regions, demonstrating that robust mitigation 
Strategies can significantly avoid future warming. 
SSP3-7.0 presents less warming than SSP5-8.5. In con- 
trast, SSP5-8.5 reflects the steepest temperature increase 
toward the end of the century. The data suggest that tem- 
perature increases in Asia may exceed the global average, 
indicating a higher regional sensitivity to global warming. 
Asia might experience significant warming under all SSP 
scenarios. Elevated areas are particularly vulnerable to 
the impacts of climate change (Hock et al., 2019). By the 
end of the century, the SSP1-2.6 scenario predicts moder- 
ate warming, whereas SSP3-7.0 and SSP5-8.5 project 
more severe temperature increases in the Alps. 


3.2 | Tundra climate zone in Asia 

Our analysis of the ET (Tundra) climate zone in 
Asia reveals no significant reduction in this zone during 
the early part of the century (Figure 2). However, sub- 
stantial portions of ET might shrink significantly under 
the SSP3-7.0 and SSP5-8.5 scenarios by mid-century, par- 
ticularly by the end of the century (more than 60%). The 
shrinkage of ET seriously threatens water availability in 
Central Asia (Didovets et al, 2024; Fallah & 
Rostami, 2024; Flanner et al., 2011). The Tibetan Plateau, 
functioning as a thermal source, is crucial in influencing 
the Asian monsoon (Immerzeel et al., 2010; Wu 
et al., 2012). Diminished snowpack and rising snow 
lines—a phenomenon associated with “elevation- 
dependent warming”’—can reduce albedo, leading to 
increased absorption of solar radiation and further ampli- 
fication of warming (Group, 2015). These changes have 
the potential to disrupt global atmospheric circulation 
patterns. 


3.3 | Tundra climate zone in Alps 

Figure 3 illustrates the potential disappearance of the ET 
climate class across all scenarios and models by the end 
of the century. The Tundra zone in the Alps could soon 
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lose half of its coverage under all scenarios. Moreover, by 
mid-century, two-thirds of the region classified as ET 
could vanish even under the sustainable SSP1-2.6 sce- 
nario, and nearly the entire zone could disappear under 
SSP3-7.0 and SSP5-8.5. As previously discussed, losing 
the Tundra climate zone in the Alps could impose severe 
climatic stress on Europe. 


3.4 | Shifts to other climate zones 

Due to its permafrost and short growing seasons, the 
Tundra is typically too cold to support tree growth. How- 
ever, aS temperatures rise, permafrost thaws and the 
growing season extends, enabling the invasion of boreal 
forests into Tundra areas. As shown in Figure 4, under 
the SSP3-7.0 and SSP5-8.5 scenarios, nearly half of the ET 
zone in Asia is projected to transition to the Dfb (Warm 
summer continental) climate zone by the end of the cen- 
tury. This shift transforms Tundra regions into climate 
zones with increased seasonal temperature variability. 
Such changes are likely to have profound ecological 
implications, including altered species distributions, 
changes in water availability, and disruptions to existing 
ecosystems. 

ET zone in Asia might convert significantly to the Dfc 
(Subarctic, Subalpine) climate classification across all sce- 
narios toward the middle and end of the century. This 
transformation, coupled with the transition to Dfb, indi- 
cates a substantial warming trend that is already affecting 
the region's climate. The conversion of Tundra areas into 
subarctic climates suggests that the ecological landscape 
of Asia will undergo considerable and early changes. 
These shifts could lead to the expansion of boreal forests 
into former Tundra regions, potentially disrupting ecosys- 
tems, altering hydrological cycles, and impacting 
biodiversity. 

Different than in Asia, in the Alps, the ET climate 
zone shifts only into Dfb and Dfc categories. The transi- 
tion from the ET climate classification to Dfb or Dfc in 
the Alps is particularly dramatic. The complex topogra- 
phy of the Alps directly influences its climate, and warm- 
ing can cause a significant upward shift in the tree line, 
leading to a reduction in ET-classified zones. Influenced 
by Atlantic and Mediterranean weather patterns, the 
Alps are susceptible to significant climate changes, 
mainly as snow cover decreases and glaciers retreat due 
to warming. Figure 5 indicates that by the end of the cen- 
tury, under SSP3-7.0 and SSP5-8.5 scenarios, the ET cli- 
mate zone in the Alps is expected to transition almost 
completely to Dfb (25% and 40% for SSP3-7.0 and 
SSP5-8.5) and Dfc (65% and 55% for SSP3-7.0 
and SSP5-8.5) classifications. The transition from ET to 
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FIGURE 3 Asin Figure 2 but for Alps. 


Dfc is evident across almost all scenarios and time slices, 
highlighting the regional warming trend. 


4 | DISCUSSION AND 
CONCLUSION 


Our paper contributes to a few studies on climatic zone 
changes in Asia and the Alps using CMIP6 model 
outputs (Beck et al., 2018; Thrasher et al., 2022). Using 


high-resolution Köppen classification, our study reveals a 
significant reduction in these zones over the Alps and 
Tibetan Plateau. 

The analysis of high-resolution Köppen classifica- 
tions elucidates ongoing and future shifts within the 
Tundra climate zones across both the Alps and Asia 
(answer to the first scientific question Q1). At the end 
of the century, there might be notable transitions to 
other climate classifications, such as Dfb and Dfc. The 
analysis shows that the extent of these shifts is closely 
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FIGURE 4 Five ISIMIP 
models' ensemble mean for the 
migration of Köppen categories 
from ET during 1981-2010 to 
other classes for different 
scenarios and time slices 

over Asia. 
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FIGURE 5 Asin Figure 4 but over Alps. 


linked to the intensity of global warming. The results 
highlight that these shifts are already underway and 
are expected to accelerate, particularly in high-altitude 
regions, where even slight temperature increases can 
lead to catastrophic outcomes. 


EE Dsb EE Dwb Dfc EE Dsc Mmm Dwc 


We conclude that surpassing a 2°C temperature 
increase from the 1981-2010 baseline will profoundly 
impact the geographical distribution and stability of the 
Tundra climate zone in the Alps and Asia (answer to 
the second scientific question, Q2). The results indicate 
that a 2°C temperature rise would likely eliminate the 
Tundra climate in the Alps and reduce its extent by over 
70% in Asia, thereby increasing the significant future 
risk, particularly in scenarios involving higher green- 
house gas emissions. 

Köppen climate classification system provides a valu- 
able framework for categorizing global climate zones. 
However, it has several caveats: Despite applying bias 
correction and statistical downscaling, projection uncer- 
tainty remains inherent in climate output. By integrating 
outputs from five structurally independent models with 
varied climate sensitivities, we embrace a range of possi- 
ble futures, thereby incorporating inherent uncertainty 
into our results. Our approach does not claim high preci- 
sion in future projections but seeks to provide fair possi- 
ble estimates given the current state of climate science. 

Increasing the number of simulations may enlarge 
the sample size of potential future climate states. Another 
approach would be to include dynamically downscaled 
models. However, there are still very few such simula- 
tions for Asia, especially the ones driven by CMIP6 
models. That could help address some model uncer- 
tainties that statistical downscaling methods may not 
capture. It is also essential to consider the limitations of 
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statistical methods, including issues related to stationar- 
ity and the fact that such approaches need high-quality 
gridded observational datasets. Such datasets might be 
available for the Alps, but Asia still has a massive gap. 
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